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Abstract 

It is a deceptively simple question to ask how acoustic disturbances propagate 

in a non-homogeneous flowing fluid. If the fluid is barotropic and inviscid, and 

the flow is irrotational (though it may have an arbitrary time dependence), 

then the equation of motion for the velocity potential describing a sound wave 

can be put in the (3 + 1) -dimensional form Atjj = {y/—g g^ v d v ip) jyj—g = 

0. The acoustic metric g^ u (t,x) governing the propagation of sound depends 

on the density, flow velocity, and local speed of sound. Even though the 

underlying fluid dynamics is Newtonian, non-relativistic, and takes place in 

flat space + time, the fluctuations (sound waves) are governed by a Lorentzian 

spacetime geometry. 

43.20.+g, 02.40.-k, 03.40.-t, 47.10.+g 
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I. INTRODUCTION 



It is well known that for a static homogeneous inviscid fluid the propagation of sound 
waves is governed by the simple equation 

dfy = c 2 V 2 rb. (1) 

(c = speed of sound.) It is a deceptively simple question to ask what happens if the fluid is 
non-homogeneous, in motion, possibly even in non-steady motion. If the fluid is barotropic 
and inviscid, and the flow is irrotational (though it may have an arbitrary time dependence) 
then I shall show that the equation of motion for the velocity potential describing an acoustic 
disturbance can be put in the (3 + l)-dimensional form 




The propagation of sound is governed by the acoustic metric g^it, x). This acoustic metric 
describes a Lorentzian (pseudo-Riemannian) geometry and depends on the density, velocity 
of flow, and local speed of sound in the fluid. Specifically 

(3) 

In general, when the fluid is non-homogeneous and flowing, the acoustic Riemann tensor 
associated with this Lorentzian metric will be nonzero. It is quite remarkable that even 
though the underlying fluid dynamics is Newtonian, nonrelativistic, and takes place in flat 
space + time, the fluctuations (sound waves) are governed by a curved Lorentzian (pseudo- 
Riemannian) geometry. This connection between fluid dynamics and techniques more com- 
monly encountered in the context of general relativity opens up many opportunities for 
cross-pollination between the two fields. 
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II. FLUID DYNAMICS 
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A. Fundamental equations 



The fundamental equations of fluid dynamics [|T]-|3|] are the equation of continuity 

d t p + V ■ (pv) = 0, (4) 

and Euler's equation 

p a ^- = p( dt v+{v-V)v)=F. (5) 

Start by assuming the fluid to be inviscid, with the only forces present being those due to 
pressure and Newtonian gravity 

F = -Vp - pV0. (6) 

Standard manipulations yield 

d t v = v x (V x v) - -JJp- V(iv 2 + <p). (7) 

Take the flow to be irrotational, introducing the velocity potential ip such that v = —Vip. 
Take the fluid to be barotropic (p is a function of p only). Then define 

fP dv' 1 

C(P) = / T~i\ ' so that V C = -Vp. (8) 
Jo p{p') p 

Euler's equation reduces to 

-d^ + c + i(WO 2 + = o. (9) 



B. Fluctuations 

Linearize the equations of motion around some assumed background (po,Po, ipo) by setting 
P = Po + e Pi, P — Po + epi, and t/j = + The gravitational potential is taken to 
be fixed and external. Sound is defined to be these linearized fluctuations in the dynamical 
quantities. The linearized continuity equation reads 
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d t pi + V • (pivo + p vi) = 0, 



(10) 



while from the Euler equation, using C = Co + e (Pi/Po), 



- d t ip\ + — - v ■ Vipi = 0. 
Po 



(11) 



Rearranging 



Pi = Po{d t ipi + vq ■ V^i). 



(12) 



Substitute this linearized Euler equation into the linearized equation of continuity. Use 
pi = [dpjdp)p\. One obtains, up to an overall sign, 

- d t ^p {dti)i + • W>i)) + V • (poVV>i - |Wo (d t ^ + v ■ W>i)) = 0. (13) 

This wave equation describes the propagation of the scalar potential tpi, and thereby com- 
pletely determines the quantities p\ and p 1 . The background fields p , p and v are per- 
mitted to have arbitrary temporal and spatial dependencies. Now, written in this form, the 
physical import of the wave equation is somewhat less than pellucid. Define 1/c 2 = dp/ dp, 
and construct the symmetric 4x4 matrix 



-1 



~v J o 



(14) 



— Vq : (c 2 5 tj — VqvI) 
Then, using four dimensional coordinates x M = (t, x l ) the wave equation is easily rewritten 

as 



o,(rdM = o. 



(15) 



This remarkably compact formulation is much more promising. The remaining steps are a 
straightforward application of the techniques of curved space Lorentzian geometry. 
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III. LORENTZIAN GEOMETRY 



In any Lorentzian (pseudo-Riemannian) manifold the curved space scalar d'Alembertian 
is given in terms of the metric g^ u (t, x) by 

1 



-A, 



-g 



-g gTdM. 



(16) 



The inverse metric, g^ u (t, x) is pointwise the matrix inverse of g^t, x), while g = det(<7 M „). 
Thus we can rewrite our physically derived wave equation in terms of the d'Alembertian 
provided we identify 



-g <r = r 



(17) 



This implies 



9 = det(D 
= (po/c 2 ) 4 [-l-(c 2 -, 2 )-, 2 ](c 2 ) 2 = -^/c 2 . 



(18) 



We can thus pick off the coefficients of the inverse metric 

-1 ; 



l 

Poc 



(19) 



— v o : (c 2 5 lj — VqVq) 
One could now determine the metric itself by inverting this 4x4 matrix. On the other hand 
it is even easier to recognize that we have in front of us an example of the Arnowitt-Deser- 
Misner split of a (3 + l)-dimensional Lorentzian spacetime metric into space + time, more 
commonly used in discussing the initial value data in general relativity (see, for example, 
pp 505-508). The metric is 



9w 



c 



5 ij 



(20) 



Observe that the signature of this metric is in fact (— , +, +, +), as it should be. 



It should be emphasised that there are two distinct metrics relevant to the current 
discussion. The physical spacetime metric is just the usual flat metric of Minkowski space 
7]fj, u = (diag[— c 2 ^, 1, 1, l])nv (coo = speed of light.) The fluid particles couple only to the 
physical metric 7]^ v . In fact the fluid motion is completely non-relativistic - ||i>o|| ^ 
Cqo. Sound waves on the other hand, do not "see" the physical metric at all. Acoustic 
perturbations couple only to the acoustic metric g^ v . The geometry determined by the 
acoustic metric does however inherit some key properties from the existence of the underlying 
flat physical metric. 

For instance, the topology of the manifold does not depend on the particular metric 
considered. The acoustic geometry inherits the underlying topology of the physical metric 
- 3? 4 with possibly a few regions excised (due to imposed boundary conditions). 

Furthermore the acoustic geometry automatically inherits the property of "stable causal- 
ity" f|. Note that ^"(V^t) (V„t) = -l/(p c) < 0. This precludes some of the more 
entertaining pathologies that sometimes arise in general relativity. 

Another concept that translates immediately is that of an "ergo-region" . Consider inte- 
gral curves of the vector (d/dty = (1,0,0,0)^. Then g ^ u (d / dty (d / dt) u = g tt = -[c 2 - ug]. 
This changes sign when | \v \ | > c. Thus any region of supersonic flow is an ergo-region. The 
analogue of this behaviour in general relativity is the ergosphere surrounding any spinning 
black hole — it is a region where space "moves" with superluminal velocity relative to the 
fixed stars. 

Observe that in a completely general Lorentzian geometry the metric has 6 degrees of 
freedom per point in spacetime. (4x4 symmetric matrix =^> 10 independent components; 
then subtract 4 coordinate conditions). In contrast, the acoustic metric is more constrained. 
Being specified completely by the three scalars ip (t,x), Po(t,x), and c(t,x), the acoustic 
metric has only 3 degrees of freedom per point in spacetime. 

A point of notation: Where the general relativist uses the word "stationary" the fluid 
dynamicist uses the phrase "steady flow" . Where the general relativist uses the word "static" 
the fluid mechanic would translate this as "fluid at rest" . 
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The analogies I am invoking between acoustics in fluids and general relativity are very 
deep and very powerful — there is a lot of mathematical machinery available for use. 



IV. GEOMETRIC ACOUSTICS 



Taking the short wavelength/high frequency limit to obtain geometrical acoustics is now 
easy. Sound rays (phonons) follow the null geodesies of the acoustic metric. Compare this to 
general relativity where in the geometrical optics approximation light rays (photons) follow 
null geodesies of the physical spacetime metric. Since null geodesies are insensitive to any 
overall conformal factor in the metric |||7j one might as well simplify life by considering the 
metric 



-Vn 



5 ij 



(21) 



Thus, in the geometric acoustics limit, sound propagation is insensitive to the density of the 
fluid, and depends only on the local speed of sound and the velocity of the fluid. It is only 
for specifically wave related properties that the density of the medium becomes important. 

One can rephrase this in a language more familiar to the acoustics community. Take 
ipx ~ ae lip . Then, neglecting variations in the amplitude a, the wave equation reduces to the 
Eikonal equation 



W v d„<p d v <p = 0. 



(22) 



This Eikonal equation is blatantly insensitive to any overall multiplicative prefactor. 

As a sanity check on the formalism, let the null geodesic be parameterized by X^(i) = 
(t,x(t)). Then the null condition implies 

dX* dX v 



h 



flU 



dt dt 



dx 







vi) - 2v. 



■ dx 1 dx 1 dx 1 



dt 



+ 



dt dt 







-dt- V ° 



c. 



(23) 



7 



Here the norm is taken in the flat physical metric. This has the obvious interpretation that 
the ray travels at the speed of sound relative to the moving medium. 

If the geometry is stationary, one can do slightly better. Let X^(s) = (t(s); x(s)) be some 
null path from x\ to X2 parameterized in terms of physical arc length (i.e. \\dx/ds\\ = 1). 
Then the condition for the path to be null (though not yet necessarily a null geodesic) is 

Solving the quadratic 



' d£ 



\ds J c 2 — v\ 

The total time taken to traverse the path is thus 



T[ 7 ] = r\dt/ds)ds 



(25) 



{ - v D ds2 + (vhdx*) 2 - vldx l }/(c 2 - vl). (26) 

Extremizing the total time taken is Fermat's principle for sound rays. One has thus checked 
the formalism for stationary geometries (steady flow) by reproducing the discussion on p 
262 of Landau and Lifshitz . 



V. DISCUSSION 
A. Limitations 

The derivation of the wave equation made two key assumptions — the flow is irrotational 
flow and the fluid is barotropic. 

The d'Alembertian equation of motion for acoustic disturbances, though derived only 
under the assumption of irrotational flow, continues to make perfectly good sense in its own 
right if the background velocity field Vq is given some vorticity. This leads one to hope that 
it might be possible to find a suitable generalization of the present derivation that might 
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work for flows with nonzero vorticity. In this regard, note that if the vorticity is everywhere 
confined to thin vortex filaments, the present derivation already works everywhere outside 
the vortex filaments themselves. 

The restriction to a barotropic fluid (p a function of p only) is in fact related to issues 
of vorticity. Examples of barotropic fluids are: 

• Isothermal fluids subject to isothermal perturbations. 

• Fluids in convective equilibrium subject to adiabatic perturbations. 

See for example fLj, §311, pp 547-548, and §313 pp 554-556. Failure of the barotropic condi- 
tion implies that the perturbations cannot be vorticity free and requires more sophisticated 
analysis. 

B. Precursors 

It is perhaps surprising that anything new can be said about so venerable a subject 
as fluid dynamics. Certainly there are precursors to the discussion of this letter in the 
fluid dynamics literature. For instance, take the background to be static, so that vq = 0, 
while d t po = = d t p , though p and hence c are permitted to retain arbitrary spatial 
dependencies. Then the wave equation derived in this letter reduces to 

«9^ = c 2 -V-(p o VV0. (27) 
Po 

This equation is in fact well known. It is equivalent, for instance to eq. (13) of §313 of 
Lamb's classic Hydrodynamics [[J. See also eq. (1.4.5) of the recent book by DeSanto |§. 
The superficially similar wave equations discussed by Landau and Lifshitz [Q] (see §74, 
eq. (74.1)), and by Skudrzyk (see p 282), utilize different physical assumptions concerning 
the behaviour of the fluid. The novelty I am claiming for the current letter is firstly, inclusion 
of nonzero background velocities and arbitrary time dependencies, and secondly and more 
importantly, the interpretation of these results in terms of Lorentzian geometry. 
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VI. CONCLUSIONS 



I have shown that acoustic waves in an inviscid fluid can (under the assumptions of irro- 
tational barotropic flow) be described by the scalar d'Alembertian of a suitable Lorentzian 
geometry. For inhomogeneous flows this Lorentzian geometry will exhibit nonzero curvature. 

Prior to this observation, Lorentzian geometries have been of interest to physics only 
within the confines of Einstein's theory of gravitation (general relativity). A large quantity 
of technical mathematical machinery currently utilized only within the context of general 
relativity may thus become of interest to the fluid dynamics community. 

On the other hand, the results of this letter give the general relativists a very down to 
earth physical model for certain classes of Lorentzian geometry. 

Particularly intriguing is the fact that while the underlying physics of fluid dynamics is 
completely nonrelativistic, Newtonian, and sharply separates the notions of space and time, 
one nevertheless sees that the fluctuations couple to a full-fledged Lorentzian spacetime. 

Note Added: After this paper was submitted for publication I was informed that similar 



results can be found in the interesting but little-known work of Unruh fllOfl . In that Let- 
ter, Unruh investigated the acoustic equivalent of Hawking radiation arising from the fluid 
dynamical analogue of a black hole. I wish to thank Ted Jacobson and John Friedman for 
bringing this reference to my attention. Further work on acoustic black holes may be found 
in gl|. 
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